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Abstract 

We present a detailed discussion of the collinear subtraction terms needed to estab- 
lish a massive variable-flavour-number scheme for the one-particle inclusive produc- 
tion of heavy quarks in hadronic collisions. The subtraction terms are computed 
by convoluting appropriate partonic cross sections with perturbative parton distri- 
bution and fragmentation functions relying on the method of mass factorization. 
We find (with one minor exception) complete agreement with the subtraction terms 
obtained in a previous publication by comparing the zero-mass limit of a fixed-order 
calculation with the genuine massles results in the MS scheme. This presentation 
will be useful for extending the massive variable-flavour-number scheme to other 
processes. 
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1 Introduction 



Heavy-quark production in highly energetic e + e~, 77, jp, ep and pp collisions has at- 
tracted much interest in the past few years, both experimentally and theoretically. Heavy 
quarks are those with masses m ^> Aqcd so that a s (m) <C 1, where a s (fiR) is the strong- 
coupling constant at renormalization scale //r. According to this definition, the charm, 
bottom and top quarks (c, b, t) are heavy whereas the up, down and strange quarks (u, d, s) 
are light. Since a s (m) <C 1, heavy-quark production is a calculable process in perturba- 
tive QCD. The heavy-quark mass acts as a cutoff for initial- and final-state collinear 
singularities and sets the scale for the perturbative expansion in a s . 

On this basis, most of the next-to-leading-order (NLO) QCD calculations have been per- 
formed in the past pQEUHlll]. Corresponding results are reliable when m is the only large 
scale, as for example in calculations of the total cross section or if any additional scale, 
for example the transverse momentum px of the produced heavy quark in 77, jp and pp 
reactions or the lepton momentum transfer Q in deep-inelastic ep scattering (DIS), is not 
much larger than the mass m. However, when p? (or Q) is much larger than the mass, 
large logarithms \n{p\,/m 2 ) or ln(Q 2 /m 2 ) arise to all orders, so that fixed-order perturba- 
tion theory is no longer valid. As is well known, these logarithms can be resummed and, 
this way, the perturbation series can be improved. 

The isolation and resummation of large logarithms is similar to the conventional mass- 
less parton model approach, where initial- or final-state collinear singularities are ab- 
sorbed into parton distribution functions (PDF) of the incoming hadrons or photons and 
into fragmentation functions (FF) for the produced light hadrons (or photons), respec- 
tively. Therefore, this approach is usually referred to as zero-mass variable-flavour-number 
scheme (ZM-VFNS). The notion "variable flavour number" is used since, in the parton 
model, the number of active quark flavours is increased by one unit, n/ — ► n/ + 1, when 
the factorization scale crosses certain transition scales (which are usually taken to be of 
the order of the heavy-quark mass) 1 . In contrast, the fixed-order treatment, where m is 
kept as a large scale, is called the fixed-flavour-number scheme (FFNS), since the number 
of flavours in the initial state is fixed to rif — 3 (4) for charm (bottom) production. One 
can combine cross sections calculated in the FFNS after certain modifications with heavy- 
quark FFs and PDFs which contain the resummed large logarithms. This prescription is 
called the massive or general-mass VFNS (GM-VFNS) 2 . 

One might expect that the partonic cross sections calculated in the FFNS approach the 
corresponding ZM-VFNS cross sections in the limit m — > if the collinear singular terms 
proportional to ln(m 2 /s) are subtracted, i.e., the subtracted FFNS cross sections differ 
from the ZM-VFNS cross sections only by terms ~ m 2 jp\. If this was true, the FFNS 
and ZM-VFNS results for the cross sections would approach each other for p\ ^> m 2 . 
This expectation, however, is not true, as was first demonstrated by Mele and Nason 



^or a detailed discussion see the appendix in Ref. L5: and references given there. 
2 For details see, e.g., Refs. 00. 
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[S] for inclusive heavy-quark production in e + e~ annihilation at NLO (e + e~ — > QQg, 
where Q is the heavy quark). They found that the limit m — > of the cross section for 
e + e~ —> QQg and the cross section calculated with m = from the beginning (in the 
MS scheme) differ by finite terms of 0(a s ). The reason for the occurrence of these finite 
terms is the different definition of the collinear singular terms in the two approaches. In 
the ZM-VFNS calculation, the heavy-quark mass is set to zero from the beginning and 
the collinearly divergent terms are defined with the help of dimensional regularization. 
This fixes the finite terms in a specific way (in a given factorization scheme), and their 
form is inherent to the chosen regularization procedure. If, on the other hand, one starts 
with m ^ and performs the limit m — > afterwards, the finite terms can be different. In 
Ref. [S], it was shown that these additional finite terms emerging in the limit m — > out 
of the theory with m/0 can be generated in the theory for m = with MS factorization 
by convoluting this cross section with a partonic fragmentation function cIq_,q(x , /i) for 
the transition from massless to massive heavy quarks Q (the explicit form of cIq-*q(x, y) 
will be given later). 

If this interpretation of the finite terms in the theory with m ^ as partonic FF is 
generally true, then cIq^q(x, y) should be process independent and could be used in any 
other heavy-quark production process. The universality of the partonic FF has been 
confirmed by performing the same calculation as in Ref. [S] for the process 7*Q — > Qg 
[3 E3 j where 7* is a space-like virtual photon, 77 — > QQg [TT] and gg — > QQg [12] and 
showing that the finite terms are obtained from a convolution of the corresponding LO 
cross sections with cIq->q(x, /j). The process-independence of cIq->q(x, fi) was established 
on more general grounds in Ref. [HJj. Moreover, process-independent derivations of the 
partonic FFs have been performed by Ma [T3j and recently by Melnikov and Mitov [131 115]. 
who have computed the partonic FFs to 0(a 2 s ). 

The fact that the theory with m 7^ and the ZM-VFNS are related by the convolution 
of the ZM-VFNS cross section with partonic FFs has been used in several ways. In 
Ref. [Bj, cIq->q(x,Ho) was used as the initial condition, at fi = /io = 0(m), for the 
calculation of cIq->q(x, (jl) at an arbitrary scale \i with the standard evolution equation. 
Later, Cacciari and Greco calculated with the same procedure the cross section for heavy- 
quark production in pp and pp collisions from the ZM-VFNS cross section supplemented 
with evolved (Iq^q^x, fi) as a function of pt [II]. Partonic FFs used together with a zero- 
mass hard-scattering cross section have subsequently been applied also to heavy-quark 
production in 77 [TH] and jp [E3 ED] processes. In Ref. [20] , the approach was generalized 
to the reaction 7 + p — > D* + X. The transition c — > D* was described by a FF containing 
besides a non-perturbative contribution the purely perturbative partonic FF. The non- 
perturbative part was described by a function containing two parameters which were fixed 
by comparison to experimental data for e + + e~ — > D* + X. In Refs. [T7] ITH] El I2H], the 
perturbative FF approach was motivated by the requirement to match the ZM-VFNS as 
close as possible to the m/0 theory. This could be achieved since at small pt = 0(m) 
the evolution of cIq_+q(x, /x ) was not yet effective and, therefore, dq^q was just taking 
account of the difference of the two theories. However, terms proportional to m 2 jp\ are 
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not included in this way. 

The so-called FONLL approach [2U E21 ESI EH ES] is an attempt to repair this deficiency. 
There, the ZM-VFNS with perturbative FFs together with a non-perturbative component 
was combined with the FFNS with rif = 3 (4) for charm (bottom) production, introducing 
a pt dependent suppression factor by hand. In addition, rr? jp\ terms have been included 
in extensions of the ACOT scheme [2EI EZj to one-particle inclusive production of D 
mesons in charged- current and neutral-current DIS [TO] EH]- In Ref. [29 , the ACOT 
scheme has been applied to one-particle inclusive heavy-quark production in pp collisions. 

Instead of incorporating the finite terms cIq-+q(x, [i) into the initial conditions of the 
perturbative FFs at fi = fiQ = 0(m), one can take this difference also into account by 
switching to a new factorization scheme, which we call the massive factorization scheme. 
In this scheme, starting from the ZM-VFNS, one adjusts the factorization of the final- 
state collinear singularities associated with the massive quarks in such a way that it 
matches the massive calculation in the limit m — > 0. Of course, the hard-scattering 
cross sections of any other process for inclusive D* production must be transformed to 
the new scheme as well. This is particularly important for the reaction e + + e~ — > 
D* + X from which the information on the non-perturbative FF for c — > D* is obtained 
by comparison to experimental data. So far, calculations in this massive factorization 
scheme were performed for 7 + p — > D* + X in Ref . [SO] , where also fits of the new FFs for 
c — > D* have been presented (for a comparison of FFs in the massive and the MS scheme, 
see Ref. [SI]). 

The simplest way to connect the truly massless cross sections in the MS scheme with 
the massive cross sections is to subtract the finite pieces cIq^q(x, /io) from the massive 
theory. In this way, one can incorporate also the m 2 jp\ terms, as given in the massive 
theory, with the advantage that the massive theory approaches the ZM-VFNS theory in 
the limit pr —> 00 or m — > 0. In addition, by including also the terms proportional to 
In m 2 contained in cIq-^q(x , /i) one can obtain not only the finite subtraction terms but 
also the terms needed for a transition to a new factorization scale. This approach has 
been applied to-f + -f^D* + X [TTJ E2], to 7 + p -> D* + X [33] and top + p^ D* + X 
[T2] . In particular in Ref. ^2] , we obtained the finite subtraction terms by comparing the 
cross sections of the massive theory, worked out by Bojak and Stratmann [SUES], in the 
limit m — > with the cross sections in the genuine massless theory in the MS factorization 
scheme as deduced by Aversa et al. jSH] in a form which is equivalent to the convolution 
of the massless cross section with cIq-,q(x, fi). 

We are going to present details of this quite involved calculation in this paper. The 
purpose is, on the one hand, to exactly demonstrate that all the subtraction terms are 
generated by the convolution with partonic FFs, at NLO just with cIq^q(x, (/,). On the 
other hand, we hope that the detailed presentation will show how the calculation carries 
over to other processes a + b — > D* + X . Since heavy-quark production in hadron-hadron 
collisions is the most complex case, we shall concentrate on this particular process. Some 
results will also be directly relevant for heavy-quark production in 77 and ^p processes. 
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Figure 1: Feynman diagrams for the LO gluon-gluon fusion process g + g — > Q + Q. 



Figure 2: The LO quark-antiquark annihilation process q + q — > Q + Q. 

The outline of this paper is as follows. In Sec. 121 we consider heavy-quark production in 
hadronic collisions, introduce the notation and review the derivation of the subtraction 
terms in Ref. ^2]- In Sec. El we present the convolution formulas, from which, in Sec. 

01 the various subtraction terms are calculated and compared with the results in Ref. 
|12j . Section El contains a summary and some concluding remarks. The subprocess cross 
sections needed in the convolutions have been collected in App.[X]for convenience. 

2 Hadroproduction of heavy quarks 

In the FFNS, the following partonic subprocesses contribute to p + p — > H + X in leading 
order (LO) and NLO, where H = D, D*, B ... is a heavy meson: 

1. g(ki) +g(k 2 ) — > Q(pi) +Q(p 2 ) + [&(P3)] ? where Q = c,b denotes a heavy quark. The 
LO Feynman diagrams are shown in Fig. ^ 

2. q(ki) + q(k 2 ) — > Q(pi) + Q(p 2 ) + [s'Cps)]- in LO, there is one Feynman diagram, 
which is shown in Fig. 121 



3. g(ki) + q(k 2 ) -> Q( Pl ) + Q{p 2 ) + q(p 3 ) and g^h) + q(k 2 ) -> Qfa) + Q{p 2 ) + q(p 3 ) 



contribute at NLO. The Feynman diagrams for these processes, as well as those for 
the NLO contributions of gg — > QQg and qq — ► QQg, can be found in App. El 




5 



Our aim is to calculate differential cross sections with an observed heavy quark Q of 
momentum p\. Therefore we define the following invariants 

s = (k 1 + k 2 ) 2 , 

ti — t — m 2 = {k\ — pi) 2 — m 2 , 

U\ = u — m 2 = [k 2 — pi) 2 — m 2 , (1) 

and 

s 2 = S 2 - m 2 = (ki + k 2 - pi) 2 -m 2 = s + ti + u\ (2) 

with s + t\ + U\ = at LO, where p 3 = 0. As usual, we introduce the dimensionless 
variables v and w by 

v = l + —, w = — , (3) 

S S + t\ 

so that 

t\ = — s(l — v), ui = —svw . (4) 

In LO, we have w — 1. 

In a recent publication [12], we have presented a NLO calculation for the inclusive produc- 
tion of D* mesons in pp collisions including heavy-quark mass effects in the hard-scattering 
cross sections. The following procedure has been adopted [HI [22] ( see a ^ so R- e fe- |71 137]): 

(i) We have derived the zero-mass limit of the cross sections in the massive FFNS with 
rif — 3 EH] for the partonic subprocesses listed above only keeping m as a 
regulator in logarithms In {m 2 / s). Special care was required in order to recover the 
distributions 5(1— w), (1/(1 — w)) + and (ln(l —w)/(l — w)) + occurring in the mass- 
less MS calculation, see, e.g., Eq. (12) in Ref. [12]. The result, generically denoted 
linim^o da(m), contains mass singular logarithms ln(m 2 ), but collinear singularities 
associated with light quarks and gluons are already subtracted in da(m). 

(ii) Then we have compared the massless limit with the corresponding hard-scattering 
cross sections in the genuine massless MS calculation in order to identify appropriate 
subtraction terms. Generically, one can write 

da sub = lim da(m) - da m , (5) 

where da^jg is a hard-scattering cross section in the genuine MS calculation. 

(iii) The desired massive hard-scattering cross sections have then been constructed by 
removing the subtraction terms from the massive cross sections in the fixed-order 
theory: 

do-(m) = da(m) -da sub . (6) 

By this procedure, the collinear logarithms ln(m 2 /s) along with finite terms which 
are independent of the heavy-quark mass are subtracted from da(m). On the other 
hand, all finite mass terms of the form (m 2 /p\) n (with an integer n) are kept in the 
hard-scattering cross sections. 
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(iv) Contributions with charm quarks in the initial state have been included in the mass- 
less approach. It can be shown that neglecting the corresponding heavy-quark mass 
terms corresponds to a convenient choice of scheme (S-ACOT scheme [38J) which 
does not imply any loss of precision. In fact, the error which is made is of the same 
order as the error of the factorization formula, as has been discussed in the con- 
text of heavy-quark production in deep inelastic scattering [3HJ EH] • Obviously, this 
rule is of great practical importance since the existing massless results for the hard- 
scattering cross sections can simply be used, whereas their massive analogues 
are unknown and would require a dedicated calculation of these processes. 3 

Note that also the FONLL calculation in Ref. [21] has been constructed with the help of 
the zero-mass limit of the fixed-order calculation in Refs. QjJ E] . On the other hand, in the 
GM-VFNS of Ref. [2H], the collinear subtractions have been obtained using the methods 
of mass factorization in a massive regularization scheme. In this approach, the subtraction 
terms are computed by convolutions of appropriate subprocesses with universal partonic 
PDFs and FFs. However, the discussion in Ref. [211] is rather generic without presenting 
many details. It is the purpose of this paper to provide a detailed description of the 
derivation of the collinear subtraction terms using the convolution method and to compare 
with the results obtained in our previous publication [T2] . 



3 Mass factorization with massive quarks 

The starting point in our approach is the basic factorization formula at the partonic level: 

da(a + b^Q + X) = fa^Xi) f b ^j(x 2 ) ® da(i +j^k + X)® d k _> Q (z) , (7) 

where da denote partonic cross sections (with singularities due to light-quark and gluon 
lines already subtracted via conventional mass factorization [30]) and da are IR-safe hard- 
scattering cross sections which are free of logarithms of the heavy-quark mass. The indices 
a, b, and i, j, k denote partons, and a sum over double indices is implied here and in 
the following. All logarithms of the heavy-quark mass (i.e. the mass singularities in 
the zero-mass limit) are contained in the partonic distribution functions / a _>j and in the 
partonic fragmentation functions dk^Q- The convolution ® denotes the usual convolution 
integral and will be specified below in Eqs. (JHJ), and (|HU|) . Equation (J2J) reflects the 
fact that the partonic cross sections da can be factorized into process-dependent IR-safe 
hard-scattering cross sections da, which are well-behaved and finite in the limit m — > 0, 
and universal (process-independent) partonic PDFs f a _i and partonic FFs d^q. 

3 For deep inelastic scattering, massive-quark-initiated coefficients have been obtained in Refs. 
the results for this simpler case are already quite involved. 
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Equation (JJJ) can be expanded in powers of a s with the help of 

fo^iM = 5 ia 5(l - Xl ) + + /«. + ..., 

h^{x 2 ) = s jb s(i - x 2 ) + fi% + /£. + . . . , 

dfc-oCr) = - z) + 4l Q + 41 Q + . . . , (8) 

da = d^°) + d<r« + d<7< 2 > + . . . , 
da = d*<°> + da« + da® + . . . . 

For the partonic PDFs and FFs, the superscript denotes the order of a s . For the cross 
sections, it indicates the relative order in a s with respect to the Born cross sections. The 
expansion of Eq. (JZj) can be used to determine order by order the relation between the 
hard-scattering and partonic cross sections. Up to NLO, one finds 

da (0) (a + b -> Q + X) = da (0) (a + b -» Q + X) = da (0) (a + b -> Q + X) , (9) 
da (1) (a + b -> Q + X) = da (1) (a + 6 -» Q + X) 

-/Si(^i) ® da (0) (« + b -> Q + X) 

-A%(x 2 ) ® da(°)(a + j - Q + X) (10) 

-da(°)(a + b -> fc + X) ® 4i«(«) • 

The three convolutions in Eq. ()10|) can be identified with the subtraction term d<r sub in 
Eq. ©. 

The factorization in Eq. (JJJ) has to be defined at a definite energy or momentum scale 
which enters as an argument into the PDFs, FFs and da. We denote the factorization 
scales by /ip for initial-state factorization (entering the PDFs) and by \J F for final-state 
factorization (entering the FFs). The renormalization scale will be called 



3.1 Partonic parton distribution and fragmentation functions 



The functions for the initial state are given in the MS scheme 4 , keeping the heavy- 



quark mass as a regulator for the collinear divergences, by 



2tt 



(i) 



[X, fMR, IJL F ) = ^ C F 



l+x 2 



X 



/4 



ln^- -21n(l-a;) - 1 



'111 



2tt 3 /i "m 2 



Note that it is assumed that the MS scheme is defined in the conventional way where photons and 
gluons have d — 2 degrees of freedom (where d is the number of space-time dimensions). Furthermore, 
subtractions /y ® da^ are performed with subprocess cross sections calculated in d dimensions. 
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where Pg-l q (x) = \[x 2 + (l — x) 2 ] and P q d q (x) = Cp[(l+x 2 )/(l — x)] + (appearing in fn^n) 
are the conventional (space-like) one- loop splitting functions and Tf = 1/2. The function 
fql^q(x, Hr, Hf) will not be used in the following, since heavy quarks in the initial state 
are treated as massless quarks as explained in Sec. |2] It would be present in cases where 
massive heavy quarks Q appear in the initial state as for example in Refs. [HI EH EE1 E2j • 

The functions d^X • for the final state read [HJ El EE El 



dqX Q {z, flR, fXp) 



2tt 
2tt 



12 

2 ' 



1 + z 2 



In 



m 2 



21n(l 



12) 



Generally, the splitting functions entering the partonic FFs are time-like splitting func- 
tions which are, however, identical to the space-like splitting functions at the one-loop 
level. It should be noted that the function fq\q{x, lir, Lip) in Eq. (fTTjl is of the same form 

as dg > _ > Q(x, [Xr, /x' f ) at 0(ol\) [3 EH- This will not be true at higher orders since the NLO 

space- and time-like splitting functions Pq^l q (x) are different. All the other partonic PDFs 
and FFs not listed here are zero at 0(a]). Furthermore, analogous results for processes 
involving photon splittings can be found by obvious replacements (g — *■ 7, a s — > a and 
appropriate modifications of colour factors) in Eqs. (JTTJ) and (|T2j). 

The partonic PDFs and FFs are known to order 0(a 2 ). They would be needed, together 
with the three-loop beta function of QCD, for computing subtraction terms at next-to- 
NLO (NNLO). For the initial state, the partonic PDFs at order 0(a 2 ) can be found in 
Ref. [41] (with the exception of /q_^q(^), which is unknown). Recently, also the 0(a 2 ) 
contributions to the perturbative FFs have been derived [THJ EE]- It should be noted 
that, at 0(a 2 ), all of the perturbative PDFs and FFs are non-vanishing at fip = m and 
fx' F = m, respectively. In fact, the parts proportional to logarithms of the factorization 
scale follow from the evolution equations, so that the new information obtained from the 
0(a 2 ) calculation is contained in the non-vanishing pieces at fxp,fx' F = m. 

In Sec. |H we will need the distribution dq^iz) with z — 1 — v + vw as a distribution in 
the kinematic variables v and w. This form of d!q_^Q(z) can be written as: 



d { ql Q (z) = A(v) 6(1 - w) + B(v) + C(v) + + D ^ w ) » ( 13 ) 
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(a) 



Pi fci 




Pi fci 



Li 



\-l—P\ 



(c) 



Figure 3: Sketch of kinematics of mass factorization for (a) upper incoming line (b) lower 
incoming line and (c) outgoing line. 



with 



A(v) 
B(v) 
Civ) 
D(v, w) 



C F 
C F 



-c F 
-c F 



2ll V 



/ 2 



In (3 + 4 In v) + 4(1 - In v - hfv) 



/ 2 

ln^r- - l-21ni; 
m 2 



2tt v 



2tt 



[2-v(l-w)\ 



In 



/ 2 

m 2 



1 -21nv-21n(l - w) 



3.2 Subtraction terms at NLO 

We distinguish mass factorization in the initial state and in the final state. For one-particle 
inclusive production, where one final-state particle has a fixed momentum (above, we had 
chosen pi), we have to distinguish further two cases with initial-state singularities corre- 
sponding to t- and li-channel scattering. A graphical representation of the subtraction 
terms in form of cut diagrams for all cases is shown in Fig. El These diagrams can be 
found by applying all possible cuts to internal lines of the Feynman diagrams (see App. 
|B| . The cuts are allowed if the 2 — > 2 subprocesses are kinematically possible and the 
1 — > 2 process involves the splitting into a heavy-quark line. In an axial gauge, the cut 
diagrams correspond to actual Feynman diagrams. 

3.2.1 Initial-state factorization 

In the first case with w-channel scattering (see Fig. Eta)), the collinear subtractions are 
given by 



da suh (a + b^Q + X) 



dxt fSliixuHtoUp) d<r (0) (i(xifci) +&(*&) -» Q(pi) + X) 



(14) 
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Here a + b — > Q + X stands for the one-particle inclusive partonic subprocesses (g + g — *> 
Q + X, g + g^Q + X, g + g-^Q + X, g + g^Q + X), / a _n(xi, /i^) describes the 
collinear splitting of parton 'a' into parton and i + b — > Q + X are the corresponding 
2^2 subprocesses with momenta xifci, k 2 and pi. A sum over i is implied, i.e., all 
possible splittings and subprocesses have to be taken into account. 

We define the following invariants for the subprocess: 

s = (xi^i + k 2 ) 2 = x v s , 

i\ = (x\ki — pi) 2 — m 2 = x\ti , (15) 
Mi = (k 2 - pi) 2 -m 2 = ui , 

and 

v = 1 + — = v , u> = — = — , (lb) 

S S+tt Xi 

t\ = —s{\ — v) , Ui = — svw = — sv . (17) 

For the calculation of d 2 cr sub /(dt>dw) in Eq. (JUJ), it is convenient to write the subprocess 
cross section as 

d ^(°) _ jd 2 a® _ jd *(°> 
di;du> d-Odw dt> 

with 

5(1 — u)) = wS(xi — X\) , Xi = w . (19) 

The 5-function imposes the 2-^2 process kinematics s + t\ + u\ = 0, i.e. zt> = 1, and 
implies x\ = w. The Jacobian reads 

J=§^ = i. (20) 
o(v,w) X\ 

Combining these results we find 

d 2 a® da® 



dvdw dv 



Sixt-xt), (21) 

so that the subtraction terms for initial-state factorization of the upper incoming line can 
be calculated as 



d 2 a su ^ da® 

-^^{a + b^Q + X) = fSii(x U fi R , fi F ) -^-{i + b^Q + X) 



• (22) 



In the second case with t-channel scattering (see Fig. Efb)), the collinear subtractions are 
given by 

da sub (a + b^Q + X) = j dx 2 f^{x 2 , fi R , fi F ) da® (a^h) + j (x 2 k 2 ) Q{ Pl ) + X) 

= fi%{x 2 ) ® da® (a + j^Q + X). (23) 
11 



The invariants for the subprocess are now given by 



s — (k 1 + x 2 k 2 ) 2 = x 2 s , 

i 1 = (k 1 -p 1 ) 2 -m 2 =t 1} (24) 
Wi = (x 2 k 2 - pi) 2 - m 2 = x 2 u\ , 



and 



x 2 -l + v A X 2 VW . 

v = , w = — — . 25 

x 2 x 2 -l + v 

Again, we write the subprocess cross section as in Eq. (jl8)l with 

S(l-u,) = x*™S( X2 -x 2 ), x 2 = -i^iL, (26) 
1 — v 1 — vw 

and 

o{v, w) x 2 — 1 + V 

For the 2-^2 subprocess kinematics, we have w = 1, x 2 = x 2 , v = vw and J = l/(x 2 w). 
Combining these results, we find 



d 2 a(°) v da(°) 



dvdw 1 — vw dv 



5(x 2 - x 2 ) , (28) 

s— >&2S, ii—*vw 



so that the subtraction terms for initial-state factorization of the lower incoming line can 
be calculated as 



dvdw 



{a + b^Q + X) = f£L(x 2 , fi R , (i F ) ; ^-(« + j ^ Q + X) 

J 1 — vw dv 



(29) 

V— *VV1, S— >X2S 



3.2.2 Final-state factorization 



The case shown in Fig. E(c) corresponds to factorization of singularities in the final state. 
Here the collinear subtractions are given by 

da sub (a + b -> Q + X) = [ dz d£<°> (a(h) + b{k 2 ) -> k{z~ l Pl ) + X) d%i Q (z, (i R , (i' F ) 

Jo 

= da^(a + b^k + X)®d { ^l Q (z). (30) 
The invariants for the subprocess can be defined as follows 

s = (ki + k 2 ) 2 = s , 

i 1 = (k 1 -z- 1 pi) 2 -m 2 = -t 1 , (31) 

ui = (k 2 - z^pi) 2 - m 2 = -u\ , 

z 
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and 



z-l + v „ vw 

w = . (32) 



z z — 1 + V 



As before, we write the subprocess cross section as in Eq. (|T5j) with 

5(1 — w) = vw5(z — z) , z — 1 —v + vw , (33) 

and 

J= ^ = i__v (34) 

o(v } w) z z — 1 + V 

From w — 1 one finds z — 1 — v + vw, v = vw/z and J = l/(zw). Combining these 
results, we find 

^ ,(,-,), (35) 



dvdw z dv 

so that the subtraction terms for final-state factorization can be calculated as 



—— a + 6 -► Q + X )= d\ l i Q {z, fi R , ti' F )- -rr- (a + b -> fc + X 
dwdw v z du 



(36) 



v^>vw/z, 



3.3 Scheme dependence and implementation freedom 

Before we turn to a discussion of our results for the collinear subtraction terms calculated 
according to Eqs. (J22J), fj2U|) and ((HE]) , we add some additional remarks. 

The partonic PDFs and FFs introduced in Sec. 13. II are given in the MS factorization and 
renormalization scheme. However, in the FFNS calculations of heavy-quark production 
[H El El HI EH 135] . a modification of the MS scheme has been adopted, called MS m or 
decoupling scheme [12], where divergences due to light quarks and gluons are treated in 
the MS scheme and divergences arising from heavy-quark loops are subtracted at zero 
momentum. In order to switch from the MS m to the MS scheme the following terms have 
to be added to the partonic cross sections of the fixed-order calculations (see Sec. 3 in 
Ref. EH): 

-a.M^ln^S? (q + q^Q + X), (37) 
Sir m z 



-a s ( m )^\n^fda^ (g + g - Q + X) . (38) 

In Eqs. (|3"Tj) and (jSHJ), the parts proportional to hifi R are due to the change of a s when 
going from the MS m to the MS scheme. If a^ 11 ^ (/■*#) and «s denote the strong- 

coupling constants in the MS m and MS schemes, respectively, one can derive from the 
renormalization group equation the following relationship between the couplings: 

a?'" 1 '^ = a< n 'W) f 1 " ^^ T / ln §) + °^s) ■ (39) 



,2 
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Figure 4: Feynman diagrams representing (a) fg-i g {xi) <S> da^(gg — > QQ) and (b) 
fg l X g {x2)®da^ '(gg — > QQ). The fermion loops on the external gluon lines are heavy-quark 
loops. 

The parts proportional to In [jl 2 f can be obtained by subtracting from the cross section in 
the gluon-gluon channel the term (fg-Lg(x-\) + fg 1 i g {x2)) <8> da^(g + g — > Q + Q) (see 

Fig. EJ). Since the function fg l X g (x) in Eq. (fTT]) is proportional to Dirac's delta function, 
this amounts to a simple multiplication with the Born cross section in the gluon-gluon 
channel. This subtraction term takes into account the different treatment of heavy-quark 
loop contributions to external gluon lines in the MS and the MS m schemes. The coefficients 
in Ref. are given in the MS m scheme. Changing the results in Ref. ^2] to the MS 
scheme according to Eqs. (f3*Tj) and (jHHjl has the expected effect of replacing [3^ ^ by 
/3g ^ in the coefficients dx and d\ in Eqs. (28), (29) and (55) of Ref. [T2], so that in the 
MS scheme Adi = Ad x = in Eqs. (35), (36) and (59) of Ref. 

Even fixing the factorization scheme to the MS scheme leaves some freedom in the im- 
plementation of a massive VFNS, as has been discussed for the case of deep inelastic 
scattering in Ref. jH]. Consider for example the condition 

lim (daOm) - da sub (m)) = da m , (40) 

which might be used as an attempt to define the subtraction terms, da^g is the massless 
hard-scattering cross section in the MS scheme. However, this requirement fixes the sub- 
traction term dcx sub (/i/m, m/px) only up to terms m/pT vanishing in the limit m — > 0. 
The precise treatment of such terms proportional to m/pT is not prescribed by factoriza- 
tion. The prescription in Eq. (J3J), dcr sub = linin^o der(m) — da^jg, is minimal in the sense 
that no finite mass terms are removed from the hard part in addition to the collinear 
logarithms ln(m 2 /s). 

The same is true from the point of view of mass factorization: The factorization and 
renormalization scheme unambiguously determines the partonic PDFs and FFs. However, 
the convolution prescription leaves some freedom in the choice of the integration variable 
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and, therefore, is only unique up to terms of the order m/Q (where Q is the hard scale). 
One example is the ACOT-x prescription jB] in inclusive DIS, which guarantees the correct 
threshold behaviour of the heavy-quark-initiated contributions. Furthermore, it is possible 
to retain the mass terms in the subprocess cross sections entering the convolution formulas. 
Actually, this is done in the original ACOT scheme [2E1 Elj • 



4 Subtraction terms: results 

We now present the results for the collinear subtraction terms, calculated using Eqs. 
(Unj) and (|3l)|). The universal partonic PDFs can be found in Sec. 13.11 The required 
subprocess cross sections have been listed for completeness in App. |XJ We retain the 
heavy-quark mass terms in the subprocess cross sections entering the convolutions. In 
order to compare with our results in Ref. ^2], these mass terms have to be dropped. 

In order to facilitate the comparison with our previous results, we expand the subtraction 
cross section in the following form, 

d2 °" sub at/, x a ( 1 "\ a (MX-™) 



dvdw \ 1 — w ) , \ 1 — w 

+Ac 5 \nv + Aci ln(l - w) + Ac n , (41) 

and use the abbreviations 

X = l-vw, Y = l-v + vw, Vi = i-v (i = 1,2). (42) 

4.1 Subtraction terms for g + g^Q + Q + g 

The coefficients Aq are decomposed into an Abelian and two non-Abelian parts, following 
Ref. [US]: 

Ac, = C(s) (c|Acf d + S AC + ±Ac*) , (43) 



with 



c ^ = m^w>- (44) 



There are four different cut diagrams contributing to the total result, which we discuss in 
the following. 

d&M(gg^QQ)®d%l Q (z) 

The cut diagrams are shown in Fig. EJ The cross section d<j( ) (gg —> QQ) is proportional 
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Figure 5: Feynman diagrams representing dfr(°\gg — > QQ) <8> Oq_^q(z). 



to the function 



r(v) 



1 — v 4m 2 
+ 



1 



m 



1 — v v sv(l — v) \ sv(l — v) 
which appears in the following expressions for the Aq. They are given by 



Aci 



o \ f 2 



+ In f In 



+ 1 — In v — In v 



x 2C(s)C F r(v)[C F - 0(1 - 



Ac 2 
Ac 3 



In 



1 - 21nt> x 2C(s)C F r(w)[C F - C A v(l - v)} 



-2 x 2C(s)C f t(v) [C f - C A v(l - v)] , 



where 



Acf d 



Ar, = C(.s)(C|A C r d + ^(C^-l)Av 



2w 2(2-2w + t> 2 ) 2v 2 w 4t> 

H 1 

v\ vw v\ Y 

m 2 f r 8v(3-2v) 8{2-2v + v 2 ) 

S \ V\ vw 

m A /-8v(ll - lhv + Qv 2 ) , 8vi(2 - 2w + w 2 ^ 

"I - o f 



3w 2 w; 



«1 



V 2 W 2 



8(2 + Av - 7v 2 + Av 3 ) 8v 2 (-5 + Av)w 8v 3 w 2 



f 2 W 



1 



V 



ac 



-4u + 



5ff 



i 8w 2 f i | 4t> (3 — Qv + 4t> 2 ) ^ m 



f 2 



F 



2 / l&v 
lQv + 



+ m 4 /16i>(3 - 2v) 16(2 - 2w + w 2 ) 



vw 



s 

16v 2 w 



We find that Ac£ q 



Ac 5 and, finally, 
Acio = Ac 5 



Ac n = ±Ac 5 1 - In 



Y 
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Figure 6: Feynman diagrams representing of da^°\gg — » gg) ® ed-^oW- 



— » 




Figure 7: Feynman diagrams representing of / ( 



^1, 



da<®{Qg^>Qg). 



The latter two relations can be derived from dgl^Q^) in Eq. (fT^j) inspecting the expressions 
for B(v), C(v) and D(v,w). (Note that the parts with B{y) and C(v) contribute to Acn 
and Acio, respectively, in cases where the 1/(1 — w) + is canceled by factors (1 — to).) 



Now we turn to a comparison with the results in Ref. [12], which have been derived as 
described in Sec. EI For fjf F — m (and neglecting terms of the order 0(m 2 /s)), Eqs. (JIBJ)- 
(JHHJ) are in complete agreement with Eqs. (18) and (21)-(24) in Ref. Furthermore, 
the parts proportional to ln(fi' F /m 2 ) in Eqs. and (}4Tj) are identical to Eqs. (37) and 
(38) in Ref. As for Acn in Eq. ()53|h the part proportional to ln(/i^/m 2 ) is in agree- 
ment with Eqs. (41) and (43) for the 'qed' and 'kq' parts. The 'oq' part reproduces Eq. 
(42) in Ref. ^2] after adding the contribution given in Eq. (|54j) . 

&<TW{gg^gg)®d%l Q {z) 

The cut diagram Fig. |U] only contributes to the part of Ac n proportional to C\: 

&i(l -2v + 2v 2 ) 16(1 - 3i; + 2v 2 ) 



Ac°? 



- 48^ + 



WW' 



k; 2 (7 - Uv + 8f 2 )w 16w 3 (-l + 2v)w 2 16v 4 w 3 



8«uf _ 8(3t» - 5t> 2 + 2t> 3 ) 8(7w - 6t» 2 + 2v 3 )^ ^ fi' F 2 



Y 3 



Y 2 



Y 



(54) 



The u-channel cut in the initial state described by the diagram in Fig. contributes: 
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Figure 8: Feynman diagrams representing of f^q{x2) <8> da^pQ — ► Q<?)- 



Aeff 



Ac?? 



1 + v 2 , o. , up 

1 - 2w + 2w 2 ) In ^-V 

vw m 2 

1+v 2 2 



vr w 



l-2w + 2w 2 ) In 



m 2 



-Ac? 



,°<i 



(55) 

(56) 
(57) 



ti 1 lo(x 2 )®d&M(gQ^Qg) 



Finally, we have a contribution from the w-channel cut in the initial state described by 
the diagram in Fig. |H1 



Ac?? d 



v(-l+2v) v 2 w { 2viv 2v v (3 - 4v + 2f 2 )\ j y r 



Vi 



X 3 X 2 



(5f 



Ac?? = 



/2v 4(1 - 2t> + 2v 2 ) 4(l-4w + 2w 2 ) 4wi 4u(l - 2u) 



V\VW Z 



+ 



+ 



X 2 



+ 



X 



x 



In 



m 2 



Ac?? = -Ac?? . 



(59) 
(60) 



The sum of Eqs. <^ and fl5E|), that of Eqs. 1|5E |) and ((55 )1 . and that of Eqs. (57J) and (jHOJ) 
are identical to Eqs. (44), (45), and (46) in Ref. [T2], respectively. 



4.2 Subtraction terms for g + g^Q + Q + g 

The results for the coefficients Aq have the following colour decomposition: 

= C q (s)^f (CpAcf + CUAef ) , 



with 



C fl (a) 



2Xs 



(61) 



(62) 



There are two different cut diagrams contributing to the total result. 
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Figure 9: Feynman diagrams representing of da^(qq — > QQ) <8) <^qLq 



z . 



Figure El shows the diagram with a cut in the final state. The 2^2 process cross section 
dcr(qq — ► QQ) is proportional to the function 



r,(u) = (1 -f) 2 + w 2 + 



2m 2 



(63) 



which will occur in the results below. Furthermore, for this contribution, the Ca parts 
vanish, i.e., Ac£ a = (i — 1,2,3,5,10,11). Therefore, the final results are proportional 



to the colour factor C%: 



Ac 2 
Ac 3 



— + In v \ In + 1 — In v — In 2 v 



x2C q (s)T q (v)C 2 F . 



In 



/ 2 

m 2 



-l-21nv) x 2C q (s)r q (v)C F . 



Ac 5 = 2C g (s)C 2 
Finally, we find again 



2 x 2C q (s)r q (v)C 2 , 

2vv\ 2v 2 Vi 3v — Qv 2 + 4v 3 



Y 3 



+ 



Y 2 



Y 



m 

f —2r ( 1 



Acio = Ac 5 
Acn = ±Ac 5 ( 1 - In 



(64) 

(65) 
(66) 

(67) 

(68) 
(69) 



One can observe the same structure of the results as for da^ (gg — > QQ) <g> 



in Sec. O 



qUqW given 



Now we turn again to the comparison with the results in Ref. [T2] . For fi' F = m (and 
neglecting terms of the order 0(m 2 /s)), Eqs. (|B3 j) -(|B9 J) are identical to Eqs. (51)-(54) in 
Ref. (T2j. The parts proportional to ln(fi' F /m 2 ) in Eqs. (|64|) and (|65p coincide with Eqs. 
(60) and (61) in Ref. [12]. As for Acn in Eq. ( f fiQ j) . the part proportional to ln(fi' F /m 2 ) 
is in agreement with Eq. (62) in Ref. ^2] only after including the contribution from 
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Figure 10: Feynman diagrams representing of da^-°\qq — > (7(7) (g) crZ^z). 



da(qq —>■ gg) ® d^Q(z), which will be given in the next subsection. 

da^{qq^gg)®d^l Q {z) 

The result for the cut diagram Fig. El reads 



4w 



2v(3 - 4v + 2v 2 ) 2(1 - 2w + 2v 2 ) 4v 3 w 4v 3 w 2 

V\ w V\ V\ Y 



In 



4(2 - - 4w w 
2 



m 2 

, 4vvf + 4(3v - 5v 2 + 2v 3 ) 2(9v - 12w 2 + 4t> 3 



y 3 



y 2 



y 



In 



(70) 



x 



(71) 



Since Acn in Eq. (jHUjl only has a Cp part, Eq. (|7T|l is the only contribution to the Ca 
part and hence is in agreement with Eq. (64) in Ref. [T2*] . Furthermore, it is easy to see 
that the sum of Ac^ in Eq. (J59j) . taken from the part oc ln(^ / 'm 2 ) , and Ac^ in Eq. (fTOj) 
reproduces Eq. (63) in Ref. [12]. 



4.3 Subtraction terms for g + q^Q + Q + q 

The process — ► QQg appears for the first time at NLO. It has the colour decomposition 

Aci = C OT (s) (C F Acf + CUAcf ) , (72) 



with 



a; 



2iVs 



(73) 



There are two different cut diagrams contributing to the total result. The results for the 
process gq — > QQq are the same as can be easily seen with the help of the expressions in 
Add. IOI 



d&W(gq^gq)®d™(z) 



For the cut diagram shown Fig. [TTJ we find 
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Figure 11: Feynman diagrams representing of da^(gq — > gq) ® dg\ Q (z) 





—* 






, i 


* 

o 

D 
D 



Figure 12: Feynman diagrams representing of /^q(^i) (S> da^(Qq — > Qq). 



Ac* 



Ac- 



-2w 2 



-f 2 + 



1 - 2u + 2v 2 



+ 2v 2 w 



f if 



+ 



3u - 2v 2 



In 



/ 2 



2w 2)'- 2Y j in 

v 2 {2 -4v + 3v 2 )w 2v 3 (l - 2f)w 2 2v 4 wj 3 



+ 



+ 



2Y 



(74) 
(75) 



Equations ()74j) and (|75jl are in agreement with Eqs. (69) and (70) in Ref. [T2*] . 



The contribution of the cut diagram Fig. El is given by 

Ac* = (1 - 2 W + 2 W 2 ) In ^ , (76) 

Acg = . (77) 
Equation (|75|l is identical to Eq. (68) in Ref. 



5 Conclusions and discussion 

We have presented a detailed description of the derivation of collinear subtraction terms 
with the help of mass factorization keeping the heavy-quark mass as a regulator for 
collinear divergences. As an example, we have considered heavy-quark production in 
hadron-hadron collisions, which is the most complex case. With one minor exception (see 
below), we have reproduced all the subtraction terms derived in Ref. [12] and found nice 
agreement. For a summary of the comparison, see Tabled Apart from giving additional 
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channel 




THIS pdptjl 


Ref. jig 


99 - 


■* QQg- 


Hf = m 




(18), 


(21)-(24) 






ocln^ 


PJ), dUD 


(37), 


(38) 






t 2 

ocln^ 

m 2 


(EHMQED-part], Q [KQ-part] 


(41), 


(43) 






ocln^ 

2 

oc \n^r 


(jn31)[OQ-part] + (jni 

0+©, 


(42) 

/ A A \ 

(44), 


/ a r\ ( ah\ 

(45), (46) 


qq- 


+ QQg: 


fi' F = m 


dSl-dSHD 


(51)- 


-(54) 






ocln^ 


(jSl, fSD 


(60) 


(61) 






oc In t-^r- 


(joyj) [u ^-parcj -+- ii i u)i 


(63) 








I 2 

oc In ^2- 


(EH 






gq - 


-> QQq- 














ocln^ 
ocln^ 

m 


(Hi, (USD 


(69), 
(68) 


(70) 



Table 1: Collinear subtraction terms for the partonic subprocesses g + g — ► Q + Q + g, 
q + q — > Q + Q + g and g + g^Q + Q + gin comparison with the results of Ref. [TJ] . In 
the second column, '/i^ = m' indicates that the final state factorization scale \J F has to be 

set to m in the equations in the third column. Furthermore, oc In (oc In refers to 

/ 2 2 

those parts of the equations in the third column which are proportional to In (In ^3-). 
The third and forth columns list the equation numbers for the corresponding subtraction 
terms derived in this paper and in Ref. [12] , respectively. Equations combined by a 'plus' 
sign have to be added. 

insight and providing a consistency check of our previous results, this detailed example 
will be useful for extending the GM-VFN scheme to other processes. 

Note however, that some exceptions have been found, (i) In Eq. (25) of Ref. ^2|, we have 
found an extra contribution to the coefficient Aci in the gg — > QQ channel resulting in 
a modification Aci — > Ac% — C(s)Ca\v(1 — v). This extra piece has its origin in heavy- 
quark loop contributions to the virtual corrections to the Born process gg — > QQ in Ref. 
[3*4] and has no counterpart in the results of Sec. 14.11 However, these terms are absent 
in the massless limit of the calculation in Refs. Q]. Numerically, this modification 
of Aci turned out to be negligible, (ii) In a publication by two of us [32., subtraction 
terms for the non-Abelian part of the process 75 — > QQg have been derived by comparing 
the zero-mass limit of Ref. [43] with the massless theory of Ref. jS] , which do not meet 
the expectations of mass factorization in Sec. El The subtraction terms derived this way 
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correctly describe the transition between the two theories. Obviously, if one of the theories 
uses conventions differing from the conventional MS scheme, the results will not agree with 
the subtraction terms derived via mass factorization. Whether this is indeed the reason 
for the discrepancy, can be clarified only with the help of a new full calculation. It is 
noteworthy that also for the channel 'jq — > QQq a non-vanishing result for the coefficient 
Acn (see Eq. (78) in Ref. [32J) was found, which would have been zero employing the 
methods in Sec. H3 In this case, the difference could be traced back to an error in Ref. 

Ell- 
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A Cross sections for 2^2 subprocesses 

In this appendix, we list the cross sections for all one-particle-inclusive subprocesses, 
a + b — > c + X , needed to compute the subtraction terms in Sec. |U For brevity, part X 
of the final state is not written explicitly in the following. We begin with subprocesses 
occurring in the channel g + g — > Q + Q + g, needed in Sec. 14.11 

A.l Subprocesses in g + g^Q + Q + g 

d<x (0) 1 1 

—{gg Q) = a s 2 7T ? [C F - Nv(l - v)} r(s, v) , (78) 

d<x (0) , s 2 4JV 2 1(1- xf 

{gg ->«/) = — ^- , (79) 



2 



dv y ^ s N 2 -ls x 2 

^-{gQ - Q) = ^vt-^-t \ } 2C F [C F v" + N(l - v)}—— , (80) 
dt> iV 2 — Is v v{l — v) 

^-(Qg - Q) = 4^(9Q - Q) 

= a 2 ^-^-\ l -^2C F [C F (1 - vf + Nv] —1— , (81) 
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where 



~\ v 1 - v 
t{s, v) = H h 



4m 



1 — v v 
x = -0(1 — v) . 



1 - 



m 



sv(l — v) \ sv(l — v) J ' 



(82) 
(83) 



A. 2 Subprocesses m.q + q^Q + Q + g 



dog 

d-0 



(gg -»• Q) = a 2 s 7r 



C F 1 



[l-v) 2 + v 2 + 



2m 2 



34) 



do-W 



d-0 



' ' Y V .\Y(,. 0)](Ah ' 



iV 2 - 1 s 



(85) 



A. 3 Subprocesses 'mg + q^Q + Q + q and g + g^Q + Q + g 



dog 

d-0 

dog 
dv 



(?0 
(#9 



dog 
du 

dqW 

df) 



1 1 



2iV 2 s 
da(°) 



(2 - 2v + v 2 ) [(N 2 - l)v 2 + 2iV 2 (l - 



v 2 (l — v) 



,(86) 



dt> 



(99 -> 9) 



1K->1— u 



1 1 



2iV 2 



T (l + v 2 ) [(N 2 - l)i) 2 + 2v + (N 2 - 1)] — - — , (87) 
s v[l — v) z 



, , 2 C F \ 1+V 2 

(Qq ^ Q) = a s 7i— -- 



N s 1 -v) 



(Qq - Q) = -^-(Q9 - Q) ■ 



39) 



B Feynman diagrams 



In this appendix we list the bremsstrahlung Feynman diagrams contributing at NLO to 
the process p + p — > H + X (H denotes a heavy meson, D, D*, B . . .) in the FFNS. 
They are the basis to generate the cut diagrams Figs. IBT[T2l as described in Sec. 13.21 We 
show separately the subset of Feynman diagrams for gg — > QQg which, after replacing 
one of the incoming gluons by a photon, contribute also to heavy-quark photoproduction, 
H + X. 
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(a) 



(b) 



(c) 




(d) 



Figure 13: Feynman diagrams for the NLO gluon bremsstrahlung process 7+g — > Q + Q + 
g. Reversing the heavy-quark lines yields the remaining graphs. Diagrams obtained from 
the ones shown here by replacing the photon with a gluon contribute to g + g — > Q + Q + g. 




(a) (b) (c) 




(d) (e) (f) 




Figure 14: Additional Feynman diagrams for the NLO gluon bremsstrahlung process 
g + g — > Q + Q + Replacing the photons by gluons in Fig. HH1 and reversing the 
heavy-quark lines of part (a) yields the remaining graphs. 
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(a) 



(b) 



(c) 





(d) 



(e) 



Figure 15: Feynman diagrams for the NLO gluon bremsstrahlung process q+q — > Q+Q+g. 



(a) 




(c) 



(b) 





(d) 



Figure 16: Feynman diagrams for the NLO light-quark-initiated subprocess 7 + q — ► 
Q + Q + q: "Bethe-Heitler" graphs (a) and (b), "Compton" graphs (c) and (d). Diagrams 
obtained from the ones shown here by replacing the photon with a gluon contribute to 
g + q^Q + Q + q. 




Figure 17: Additional Feynman diagram for the NLO light-quark- initiated subprocess 
g + q — > Q + Q + q. Replacing the photons by gluons in Fig. HH1 yields the remaining 
graphs. 
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